Sampling the reciprocal Coulomb potential in finite anisotropic cells.

We present a robust strategy to numerically sample the Coulomb potential in reciprocal space for periodic Born-von Karman cells of general shape. Our approach tackles two common issues of plane-wave based implementations of Coulomb integrals under periodic boundary conditions: the treatment of the singularity at the Brillouin-zone center and discretization errors, which can cause severe convergence problems in anisotropic cells, necessary for the calculation of low-dimensional systems. We apply our strategy to the Hartree-Fock and coupled cluster (CC) theories and discuss the consequences of different sampling strategies on different theories. We show that sampling the Coulomb potential via the widely used probe-charge Ewald method is unsuitable for CC calculations in anisotropic cells. To demonstrate the applicability of our developed approach, we study two representative, low-dimensional use cases: the infinite carbon chain, for which we report the first periodic CCSD(T) potential energy surface, and a surface slab of lithium hydride, for which we demonstrate the impact of different sampling strategies for calculating surface energies. We find that our Coulomb sampling strategy serves as a vital solution, addressing the critical need for improved accuracy in plane-wave based CC calculations for low-dimensional systems.

Electronic Free Energy Surface of the Nitrogen Dimer Using First-Principles Finite Temperature Electronic Structure Methods.

We use full configuration interaction and density matrix quantum Monte Carlo methods to calculate the electronic free energy surface of the nitrogen dimer within the free-energy Born-Oppenheimer approximation. As the temperature is raised from T = 0, we find a temperature regime in which the internal energy causes bond strengthening. At these temperatures, adding in the entropy contributions is required to cause the bond to gradually weaken with increasing temperature. We predict a thermally driven dissociation for the nitrogen dimer between 22,000 to 63,200 K depending on symmetries and basis set. Inclusion of more spatial and spin symmetries reduces the temperature required. The origin of these observations is explored using the structure of the density matrix at various temperatures and bond lengths.

Electronic specific heat capacities and entropies from density matrix quantum Monte Carlo using Gaussian process regression to find gradients of noisy data.

We present a machine learning approach to calculating electronic specific heat capacities for a variety of benchmark molecular systems. Our models are based on data from density matrix quantum Monte Carlo, which is a stochastic method that can calculate the electronic energy at finite temperature. As these energies typically have noise, numerical derivatives of the energy can be challenging to find reliably. In order to circumvent this problem, we use Gaussian process regression to model the energy and use analytical derivatives to produce the specific heat capacity. From there, we also calculate the entropy by numerical integration. We compare our results to cubic splines and finite differences in a variety of molecules in which Hamiltonians can be diagonalized exactly with full configuration interaction. We finally apply this method to look at larger molecules where exact diagonalization is not possible and make comparisons with more approximate ways to calculate the specific heat capacity and entropy.

How the Exchange Energy Can Affect the Power Laws Used to Extrapolate the Coupled Cluster Correlation Energy to the Thermodynamic Limit.

Finite size error is commonly removed from coupled cluster theory calculations by N-1 extrapolations over correlation energy calculations of different system sizes (N), where the N-1 scaling comes from the total energy rather than the correlation energy. However, previous studies in the quantum Monte Carlo community suggest an exchange-energy-like power law of N-2/3 should also be present in the correlation energy when using the conventional Coulomb interaction. The rationale for this is that the total energy goes as N-1 and the exchange energy goes as N-2/3; thus, the correlation energy should be a combination of these two power laws. Further, in coupled cluster theory, these power laws are related to the low G scaling of the transition structure factor, S(G), which is a property of the coupled cluster wave function calculated from the amplitudes. We show here that data from coupled cluster doubles calculations on the uniform electron gas fit a function with a low G behavior of S(G) ∼ G. The prefactor for this linear term is derived from the exchange energy to be consistent with an N-2/3 power law at large N. Incorporating the exchange structure factor into the transition structure factor results in a combined structure factor of S(G) ∼ G2, consistent with an N-1 scaling of the exchange-correlation energy. We then look for the presence of an N-2/3 power law in the energy. To do so, we first develop a plane-wave cutoff scheme with less noise than the traditional basis set used for the uniform electron gas. Then, we collect data from a wide range of electron numbers and densities to systematically test five methods using N-1 scaling, N-2/3 scaling, or combinations of both scaling behaviors. We find that power laws that incorporate both N-1 and N-2/3 scaling perform better than either alone, especially when the prefactor for N-2/3 scaling can be found from exchange energy calculations.

Machine learning for a finite size correction in periodic coupled cluster theory calculations.

We introduce a straightforward Gaussian process regression (GPR) model for the transition structure factor of metal periodic coupled cluster singles and doubles (CCSD) calculations. This is inspired by the method introduced by Liao and Grüneis for interpolating over the transition structure factor to obtain a finite size correction for CCSD [K. Liao and A. Grüneis, J. Chem. Phys. 145, 141102 (2016)] and by our own prior work using the transition structure factor to efficiently converge CCSD for metals to the thermodynamic limit [Mihm et al., Nat. Comput. Sci. 1, 801 (2021)]. In our CCSD-FS-GPR method to correct for finite size errors, we fit the structure factor to a 1D function in the momentum transfer, G. We then integrate over this function by projecting it onto a k-point mesh to obtain comparisons with extrapolated results. Results are shown for lithium, sodium, and the uniform electron gas.

Piecewise interaction picture density matrix quantum Monte Carlo.

The density matrix quantum Monte Carlo (DMQMC) set of methods stochastically samples the exact N-body density matrix for interacting electrons at finite temperature. We introduce a simple modification to the interaction picture DMQMC (IP-DMQMC) method that overcomes the limitation of only sampling one inverse temperature point at a time, instead allowing for the sampling of a temperature range within a single calculation, thereby reducing the computational cost. At the target inverse temperature, instead of ending the simulation, we incorporate a change of picture away from the interaction picture. The resulting equations of motion have piecewise functions and use the interaction picture in the first phase of a simulation, followed by the application of the Bloch equation once the target inverse temperature is reached. We find that the performance of this method is similar to or better than the DMQMC and IP-DMQMC algorithms in a variety of molecular test systems.

The Sign Problem in Density Matrix Quantum Monte Carlo.

Density matrix quantum Monte Carlo (DMQMC) is a recently developed method for stochastically sampling the N-particle thermal density matrix to obtain exact-on-average energies for model and ab initio systems. We report a systematic numerical study of the sign problem in DMQMC based on simulations of atomic and molecular systems. In DMQMC, the density matrix is written in an outer product basis of Slater determinants. In principle, this means that DMQMC needs to sample a space that scales in the system size, N, as O[(exp(N))2]. In practice, removing the sign problem requires a total walker population that exceeds a system-dependent critical walker population (Nc), imposing limitations on both storage and compute time. We establish that Nc for DMQMC is the square of Nc for FCIQMC. By contrast, the minimum Nc in the interaction picture modification of DMQMC (IP-DMQMC) is only linearly related to the Nc for FCIQMC. We find that this difference originates from the difference in propagation of IP-DMQMC versus canonical DMQMC: the former is asymmetric, whereas the latter is symmetric. When an asymmetric mode of propagation is used in DMQMC, there is a much greater stochastic error and is thus prohibitively expensive for DMQMC without the interaction picture adaptation. Finally, we find that the equivalence between IP-DMQMC and FCIQMC seems to extend to the initiator approximation, which is often required to study larger systems with large basis sets. This suggests that IP-DMQMC offers a way to ameliorate the cost of moving between a Slater determinant space and an outer product basis.

Iron-Catalyzed Halogen Exchange of Trifluoromethyl Arenes*.

The facile production of ArCF2 X and ArCX3 from ArCF3 using catalytic iron(III)halides is reported, which constitutes the first iron-catalyzed halogen exchange for non-aromatic C-F bonds. Theoretical calculations suggest direct activation of C-F bonds by iron coordination. ArCX3 and ArCF2 X products of the reaction are synthetically valuable due to their diversification potential. In particular, chloro- and bromodifluoromethyl arenes (ArCF2 Cl, ArCF2 Br respectively) provide access to a myriad of difluoromethyl arene derivatives (ArCF2 R). To optimize for mono-halogen exchange, a statistical method called Design of Experiments was used. Optimized parameters were successfully applied to electron rich and electron deficient aromatic substrates, and to the late stage diversification of flufenoxuron, a commercial insecticide. These methods are highly practical, being run at convenient temperatures and using inexpensive common reagents.

Power Laws Used to Extrapolate the Coupled Cluster Correlation Energy to the Thermodynamic Limit.

Recent calculations using coupled cluster on solids have raised the discussion of using a N-1/3 power law to fit the correlation energy when extrapolating to the thermodynamic limit, an approach which differs from the more commonly used N-1 power law, which is, for example, often used by quantum Monte Carlo methods. In this paper, we present one way to reconcile these viewpoints. Coupled cluster doubles calculations were performed on uniform electron gases reaching system sizes of 922 electrons for an extremely wide range of densities (0.1 < rs < 100.0) to study how the correlation energy approaches the thermodynamic limit. The data were corrected for the basis set incompleteness error and use a selected twist angle approach to mitigate the finite size error from shell filling effects. Analyzing these data, we initially find that a power law of N-1/3 appears to fit the data better than a N-1 power law in the large system size limit. However, we provide an analysis of the transition structure factor showing that N-1 still applies to large system sizes and that the apparent N-1/3 power law occurs only at low N.

Accelerating convergence to the thermodynamic limit with twist angle selection applied to methods beyond many-body perturbation theory.

We recently developed a scheme to use low-cost calculations to find a single twist angle where the coupled cluster doubles energy of a single calculation matches the twist-averaged coupled cluster doubles energy in a finite unit cell. We used initiator full configuration interaction quantum Monte Carlo as an example of an exact method beyond coupled cluster doubles theory to show that this selected twist angle approach had comparable accuracy in methods beyond coupled cluster. Furthermore, at least for small system sizes, we show that the same twist angle can also be found by comparing the energy directly (at the level of second-order Moller-Plesset theory), suggesting a route toward twist angle selection, which requires minimal modification to existing codes that can perform twist averaging.

A shortcut to the thermodynamic limit for quantum many-body calculations of metals.

Computationally efficient and accurate quantum mechanical approximations to solve the many-electron Schrödinger equation are crucial for computational materials science. Methods such as coupled cluster theory show potential for widespread adoption if computational cost bottlenecks can be removed. For example, extremely dense k-point grids are required to model long-range electronic correlation effects, particularly for metals. Although these grids can be made more effective by averaging calculations over an offset (or twist angle), the resultant cost in time for coupled cluster theory is prohibitive. We show here that a single special twist angle can be found using the transition structure factor, which provides the same benefit as twist averaging with one or two orders of magnitude reduction in computational time. We demonstrate that this not only works for metal systems but also is applicable to a broader range of materials, including insulators and semiconductors.

NECI: N-Electron Configuration Interaction with an emphasis on state-of-the-art stochastic methods.

We present NECI, a state-of-the-art implementation of the Full Configuration Interaction Quantum Monte Carlo (FCIQMC) algorithm, a method based on a stochastic application of the Hamiltonian matrix on a sparse sampling of the wave function. The program utilizes a very powerful parallelization and scales efficiently to more than 24 000 central processing unit cores. In this paper, we describe the core functionalities of NECI and its recent developments. This includes the capabilities to calculate ground and excited state energies, properties via the one- and two-body reduced density matrices, as well as spectral and Green's functions for ab initio and model systems. A number of enhancements of the bare FCIQMC algorithm are available within NECI, allowing us to use a partially deterministic formulation of the algorithm, working in a spin-adapted basis or supporting transcorrelated Hamiltonians. NECI supports the FCIDUMP file format for integrals, supplying a convenient interface to numerous quantum chemistry programs, and it is licensed under GPL-3.0.

Using Density Matrix Quantum Monte Carlo for Calculating Exact-on-Average Energies for ab Initio Hamiltonians in a Finite Basis Set.

We here apply the recently developed initiator density matrix quantum Monte Carlo (i-DMQMC) to a variety of atoms and molecules in vacuum. i-DMQMC samples the exact density matrix of a Hamiltonian at finite temperature and combines the accuracy of full configuration interaction quantum Monte Carlo (FCIQMC)-full configuration interaction (FCI) or exact energies in a finite basis set-with finite temperature. In order to explore the applicability of i-DMQMC for molecular systems, we choose to study a recently developed test set by Rubenstein and co-workers: Be, H2O, and H10 at near-equilibrium and stretched geometries. We find that, for Be and H2O, i-DMQMC delivers energies with submillihartree accuracy when compared with finite temperature FCI. For H2O and both geometries of H10, we examine the difference between FT-AFQMC and i-DMQMC, which, in turn, estimates the difference in canonical versus grand canonical energies. We close with two discussions: one of simulation settings (initiator error, the interaction picture, and different basis sets), and another of energy difference calculations in the form of specific heat capacity and ionization potential calculations.

The Influence of Redox-Innocent Donor Groups in Tetradentate Ligands Derived from o-Phenylenediamine: Electronic Structure Investigations with Nickel.

The continued development of redox-active ligands requires an understanding as to how ligand modifications and related factors affect the locus of redox activity and spin density in metal complexes. Here we describe the synthesis, characterization, and electronic structure of nickel complexes containing triaryl NNNN (1) and SNNS (2) ligands derived from o-phenylenediamine. The tetradentate ligands in 1 and 2 were investigated and compared to those in metal complexes with compositionally similar ligands to determine how ligand-centered redox properties change when redox-active flanking groups are replaced with redox-innocent NMe2 or SMe. A derivative of 2 in which the phenylene backbone was replaced with ethylene (3) was also prepared to interrogate the importance of o-phenylenediamine for ligand-centered redox activity. Cyclic voltammograms collected for 1 and 2 revealed two fully reversible ligand-centered redox events. Remarkably, several quasi-reversible ligand-centered redox waves were also observed for 3 despite the absence of the o-phenylenediamine subunit. Oxidizing 1 and 2 with silver salts containing different counteranions (BF4-, OTf-, NTf2-) allowed the electrochemically generated complexes to be analyzed as a function of different oxidation states using single-crystal X-ray diffraction (XRD), EPR spectroscopy, and S K-edge X-ray absorption spectroscopy. The experimental data are corroborated by DFT calculations, and together, they reveal how the location of unpaired spin density and electronic structure in singly and doubly oxidized salts of 1 and 2 varies depending on the coordinating ability of the counteranions and exogenous ligands such as pyridine.

Fully Quantum Embedding with Density Functional Theory for Full Configuration Interaction Quantum Monte Carlo.

We here develop a fully quantum embedded version of initiator full configuration interaction quantum Monte Carlo (i-FCIQMC) and apply it to study an ionic bond (lithium hydride, LiH) and a covalent bond (hydrogen flouride, HF) physisorbed to a benzene molecule. The embedding is performed using a recently developed Huzinaga projection operator approach, which affords good synergy with i-FCIQMC by minimizing the number of orbitals in the calculation. When considering the dissociation energy of these bonds into closed-shell ionic fragments, we find that i-FCIQMC embedded in density functional theory (i-FCIQMC-in-DFT) delivers comparable accuracy with coupled cluster singles and doubles with perturbative triples embedded in density functional theory (CCSD(T)-in-DFT). In treating the bond dissociation energy curve of HF, i-FCIQMC-in-DFT has improved accuracy over CCSD(T)-in-DFT due to the presence of strong correlation. We discuss the implications of the new i-FCIQMC-in-DFT method as applied to bond breaking in catalysis.

An optimized twist angle to find the twist-averaged correlation energy applied to the uniform electron gas.

We explore an alternative to twist averaging in order to obtain more cost-effective and accurate extrapolations to the thermodynamic limit (TDL) for coupled cluster doubles (CCD) calculations. We seek a single twist angle to perform calculations at, instead of integrating over many random points or a grid. We introduce the concept of connectivity, a quantity derived from the nonzero four-index integrals in an MP2 calculation. This allows us to find a special twist angle that provides appropriate connectivity in the energy equation, which yields results comparable to full twist averaging. This special twist angle effectively makes the finite electron number CCD calculation represent the TDL more accurately, reducing the cost of twist-averaged CCD over Ns twist angles from Ns CCD calculations to Ns MP2 calculations plus one CCD calculation.

The HANDE-QMC Project: Open-Source Stochastic Quantum Chemistry from the Ground State Up.

Building on the success of Quantum Monte Carlo techniques such as diffusion Monte Carlo, alternative stochastic approaches to solve electronic structure problems have emerged over the past decade. The full configuration interaction quantum Monte Carlo (FCIQMC) method allows one to systematically approach the exact solution of such problems, for cases where very high accuracy is desired. The introduction of FCIQMC has subsequently led to the development of coupled cluster Monte Carlo (CCMC) and density matrix quantum Monte Carlo (DMQMC), allowing stochastic sampling of the coupled cluster wave function and the exact thermal density matrix, respectively. In this Article, we describe the HANDE-QMC code, an open-source implementation of FCIQMC, CCMC and DMQMC, including initiator and semistochastic adaptations. We describe our code and demonstrate its use on three example systems; a molecule (nitric oxide), a model solid (the uniform electron gas), and a real solid (diamond). An illustrative tutorial is also included.

A structural and mechanistic study of π-clamp-mediated cysteine perfluoroarylation.

Natural enzymes use local environments to tune the reactivity of amino acid side chains. In searching for small peptides with similar properties, we discovered a four-residue π-clamp motif (Phe-Cys-Pro-Phe) for regio- and chemoselective arylation of cysteine in ribosomally produced proteins. Here we report mutational, computational, and structural findings directed toward elucidating the molecular factors that drive π-clamp-mediated arylation. We show the significance of a trans conformation prolyl amide bond for the π-clamp reactivity. The π-clamp cysteine arylation reaction enthalpy of activation (ΔH‡) is significantly lower than a non-π-clamp cysteine. Solid-state NMR chemical shifts indicate the prolyl amide bond in the π-clamp motif adopts a 1:1 ratio of the cis and trans conformation, while in the reaction product Pro3 was exclusively in trans. In two structural models of the perfluoroarylated product, distinct interactions at 4.7 Å between Phe1 side chain and perfluoroaryl electrophile moiety are observed. Further, solution 19F NMR and isothermal titration calorimetry measurements suggest interactions between hydrophobic side chains in a π-clamp mutant and the perfluoroaryl probe. These studies led us to design a π-clamp mutant with an 85-fold rate enhancement. These findings will guide us toward the discovery of small reactive peptides to facilitate abiotic chemistry in water.

Speed Limit for Triplet-Exciton Transfer in Solid-State PbS Nanocrystal-Sensitized Photon Upconversion.

Hybrid interfaces combining inorganic and organic materials underpin the operation of many optoelectronic and photocatalytic systems and allow for innovative approaches to photon up- and down-conversion. However, the mechanism of exchange-mediated energy transfer of spin-triplet excitons across these interfaces remains obscure, particularly when both the macroscopic donor and acceptor are composed of many separately interacting nanoscopic moieties. Here, we study the transfer of excitons from colloidal lead sulfide (PbS) nanocrystals to the spin-triplet state of rubrene molecules. By reducing the length of the carboxylic acid ligands on the nanocrystal surface from 18 to 4 carbon atoms, thinning the effective ligand shell from 13 to 6 Å, we are able to increase the characteristic transfer rate by an order of magnitude. However, we observe that the energy transfer rate asymptotes for shorter separation distances (≤10 Å) which we attribute to the reduced Dexter coupling brought on by the increased effective dielectric constant of these solid-state devices when the aliphatic ligands are short. This implies that the shortest ligands, which hinder long-term colloidal stability, offer little advantage for energy transfer. Indeed, we find that hexanoic acid ligands are already sufficient for near-unity transfer efficiency. Using nanocrystals with these optimal-length ligands in an improved solid-state device structure, we obtain an upconversion efficiency of (7 ± 1)% with excitation at λ = 808 nm.

Salt Effect Accelerates Site-Selective Cysteine Bioconjugation.

Highly efficient and selective chemical reactions are desired. For small molecule chemistry, the reaction rate can be varied by changing the concentration, temperature, and solvent used. In contrast for large biomolecules, the reaction rate is difficult to modify by adjusting these variables because stringent biocompatible reaction conditions are required. Here we show that adding salts can change the rate constant over 4 orders of magnitude for an arylation bioconjugation reaction between a cysteine residue within a four-residue sequence (π-clamp) and a perfluoroaryl electrophile. Biocompatible ammonium sulfate significantly enhances the reaction rate without influencing the site-specificity of π-clamp mediated arylation, enabling the fast synthesis of two site-specific antibody-drug conjugates that selectively kill HER2-positive breast cancer cells. Computational and structure-reactivity studies indicate that salts may tune the reaction rate through modulating the interactions between the π-clamp hydrophobic side chains and the electrophile. On the basis of this understanding, the salt effect is extended to other bioconjugation chemistry, and a new regioselective alkylation reaction at π-clamp cysteine is developed.